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In recent years different explanations are provided for both an inflation and a recent acceleration 
in the expansion of the universe. In this Letter we show that a model of physical interest is the 
modification of general relativity with a Gauss-Bonnet term coupled to a dynamical scalar-field 
as predicted by certain versions of string theory. This construction provides a model of evolving 
dark energy that naturally explains a dynamical relaxation of the vacuum energy (gravitationally 
repulsive pressure) to a small value (exponentially close to zero) after a sufficient number of e-folds. 
The model also leads to a small deviation from the w = —1 prediction of non-evolving dark energy. 
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Inflation, or a period of accelerated expansion in the 
early universe, and cosmic acceleration at late times as 
pure gravitational dynamics are unusual within the con- 
text of general relativity (GR). Einstein back in 1917 
amended his General Theory of Relativity with a cosmo- 
logical constant A to achieve a stationary universe, but 
it was later realized that a positive A, in a cosmologi- 
cal context, gives rise to an accelerating expansion of the 
universe. In modern parlance, a positive A, more gen- 
erally, a vacuum energy, is called "dark energy" , which 
makes up about 70% of the matter-energy content of the 
present universe. The right explanation of dark energy 
is the greatest challenge faced by the current generation 
of physicists and cosmologists. 

A pure cosmological constant, A ~ [IQ^^eVY, as a 
source of dark energy does not seem much likely, as there 
is no conceivable mechanism to explain why it is so tiny 
and also why should it be comparable to the present mat- 
ter density (see the reviews 2] ) . It is very problematic 
for particle physics if A is to be interpreted as the vac- 
uum energy. If dark energy is something fundamental, 
then it is natural to attribute it to one or more scalar 
fields Q. It is not only because dark energy associated 
with the dynamics of scalar fields which are uniform in 
space provides a mechanism for generating the observed 
density perturbations 4] and a negative pressure suffi- 
cient to drive the accelerating expansion 5] , but there is 
another reason well motivated by fundamental physics. 

Gravity is attractive and thus curves spacetime whose 
dynamics is set entirely by the spacetime curvature: the 
Riemann curvature tensor. However, inflation (or cos- 
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mic acceleration) is repulsive: this is caused not by grav- 
ity altering its sense, attraction to repulsion, but due 
to an extra source, or vacuum energy, which covers ev- 
ery point in space and exerts gravitationally smooth and 
negative pressure. The fundamental scalar fields abun- 
dant in higher dimensional theories of gravity and fields, 
such as string/M theory, are such examples. Most of the 
scalars in string theory are the structure moduli associ- 
ated with the internal geometry of space, which do not di- 
rectly couple with the curvature tensor. But a scalar field 
associated with the overall size and the shape of the in- 
ternal compactification manifold generically couples with 
Riemann curvature tensor, even in four dimensions [^. 

The question that naturally arises is: what is the most 
plausible form of a four-dimensional gravitational ac- 
tion that offers a resolution to the dilemma posed by 
the current cosmic acceleration 7], or the dark energy 
problem, within a natural theoretical framework? In 
recent years, this problem has been addressed in hun- 
dreds of papers proposing various kinds of modification 
of the energy-momentum tensor in the "vacuum" (i.e., 
the source of agravitationally repulsive pressure), e.g., 
phantom field j^], f{R) gravity that adds terms propor- 
tional to inverse powers of R '(11 or the terms, or both 
these effects at once (Ifl], braneworld modifications of 
Einstein's GR ^tc. Other examples of recent inter- 
est are cosmological compactifications of string/M theory 
on a hyperbolic manifold 12] and on twisted spaces Il3| . 

A very desirable feature of a theory of scalar-tensor 
gravity is quasi-linearity: the property that the highest 
derivatives of the metric appear in the field equations 
only linearly, so as to make the theory ghost free. Inter- 
estingly, differential geometry offers a particular combi- 
nation of the curvature squared terms with such behav- 
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ior, known as the Gauss-Bonnet (GB) integrand: 
g = R^~ AR^.R^" + R^.^paR'""'". 

All versions of string theory in 10 dimensions (except 
type II) include this term as the leading order a' cor- 
rection 01 . A unique property of the string effective 
action is that the couplings are field dependent, and thus 
in principle space-time dependent. The effective action 
for the system may be taken to be 



S 



Via) + f{a)g 



(1) 



where k is the inverse Planck mass Mp^ = (SttGjv)^^^, 7 
is a coupling constant and /(cr) = X — 6^{a): the coupling 
A may be related to string coupling Qs via A ^ ^/gl- The 
numerical coefhcient 5 typically depends on the mass- 
less spectrum of every particular model IS 
phenomcnologically-motivated field potential. Our ap- 
proach is a new take on a familiar system: the Gauss- 
Bonnet combination is multiplied by a function of the 
scalar field. In four dimensions [d = 4) the GB term 
makes no contribution if /(cr) = const. It is natural 
to consider /(cr) as a dynamical variable. This follows, 
for example, from the one-loop corrected string effective 
action 6j when going from the string frame to the Ein- 
stein frame to appropriately describe the universe we ob- 
serve today, where the equivalence principle is well pre- 
served and the Newton's constant Gat is (almost) time- 
independent. A discussion on cosmological advantages 



of modified Gauss-Bonnet theory may be found in 11 
A study has been recently made in [l9| by introducing 
higher order (quartic) curvature corrections, but in a 
background of fixed modulus field a. 

The most pertinent question that we would like to ask 
is: what new features would a dynamical GB coupling 
introduce, and how can they influence a cosmological so- 
lution? In this Letter we present an exact cosmological 
solution that explains a dynamical relaxation of vacuum 
energy to a small value (exponentially close to zero) after 
a sufficient number of e- folds, leading to a small deviation 
from the w — —1 prediction of non-evolving dark energy. 

Let us consider the four-dimensional spacetime metric 
in standard Friedmann-Robertson- Walker (FRW) form: 



ds^ = -dt:^ + a{tf Y^idx'f 



(2) 



where a(t) is the scale factor of the universe. In terms of 
the following dimensionless variables 



X = (7^V2)(a/ff)^ 



y -- 

h : 



: H/H^, (3) 



the Einstein field equations obtained by varying the met- 
ric are given by 

= -3 + x + y-3{u' -2hu), (4) 
= u" + {2- h)u' -2h'u- 2(2 + h)uh 

+ 2h + 3 + x-y, (5) 

where X' = ^ = ±^ ^ a^, ^ J H dt = ln(a/ao) 
is the number of e-folds and H = - is the Hubble ex- 

a 

pansion parameter. The time evolution of a is given by 



d 
dt 



7 0-2 + 2A(o-)) = -6H (ja^) - 26, 



(6) 



where A(cr) = Via) - /(ct)^ and 6 = /(cr)^, with G = 
24i?2(ff + ff2). We will call A(cr) an effective potential. 
Eq. © may be written as x' + 2ih + 3)x + y' + 2hy — 
3ih + l)iu' ~ 2hu) = 0. 

When /((t) = (or u = 0), the equations of motion 
satisfy simple relationships: y = S + h and x — —h. From 
this we see that it is enough to introduce only one new 
parameter (e.g. the potential Via)) to explain a partic- 
ular value of slow-roll type variable h {= ad /a?' — 1) or 
the acceleration parameter a; if more than one parame- 
ter is introduced in the starting action, as such the case 
with /(a) ^ 0, then one complicates the problem. This 
remark though valid to some extend may have less sig- 
nificance for at least two reasons. Firstly, a non-trivial 
coupling between the field a and the curvature squared 
terms naturally exists in one-loop corrected string effec- 
tive action [l5|, as well as in the best motivated scalar- 
tensor theories which respect most of GR's symmetries, 
see for example Secondly, the knowledge about 

the modulus-dependent GB coupling may help one to 
construct a cosmological model where h (or the dark en- 
ergy equation of state parameter, w) is dynamical. 

While there can be a myriad of scalar field models with 
different forms of Via) and /(a), inspired by particle 
physics beyond the standard model or string theory, not 
all scalar potentials may be used to describe the universe 
that we observe today. In order for the model to work 
a scalar field should relax its potential energy after in- 
flation down to a sufficiently low value, presumably very 
close to the observed value of the dark energy in order 
to solve the cosmological constant problem. In this Let- 
ter, rather than picking up particular functional forms for 
Via) and /(cr), for example, as in Ref. (2fl|, we would like 
to obtain an exact cosmological solution, which respects 
the symmetry of the field equations. 

By solving the field equations Q-®, we can 
write X and y (and hence A(cr)) as a second or- 
der differential equation in u. Notably, K^A(cr) — 
^ [u" + i^-h)u' -2i?,h + h? + h' ^Z)u] +H'^i3 + h), 
which generally has a solution composed of homogeneous 
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and nonhomogeneous parts. We refer to the solution of 
u" + {5-h)u'- 2{8h + h'^ + h' + 3)u = 0, or equivalently, 

A{(7) = H^{3 + h) = + HH', (7) 

a homogeneous solution, which is trivially satisfied for 
f{a) = 0. Effectively, for f{a) ^ 0, we will pick V{a) = 
f{(j)Q + H'^{3 + h)/ K^, leaving us with only one arbitrary 
function, out of a:(77), ^.(77) or u{r/). 

What is the advantage of adding a new term, f{a)Q, 
if one could make an ansatz like with /(cr) = 0? An 
interesting feature of our construction is that while the 
contributions coming from both the field potential V{a) 
and the coupled GB term, f[a)G, can be large separately, 
the effective potential, A(<t), can be exponentially close 
to zero at late times, as it relaxes to a small value after a 
sufficiently large number of e- folds of expansion. More- 
over, even if we do not impose a priori the constraint 
^{'^) ~ fi'^)S = i?^(3 + h)/ we arrive at a similar ex- 
pression for various other solutions; below we will present 
one such example. 

Given a perfect fluid form for the energy momentum 
tensor of the field cr: Too ~ Pa, Tu = pcra{ty, we define 
the equation of state (EOS) parameter 

2h + 3 2q^l 
w= — = — = — - — , (8) 

Pa O 3 

where q = —ad/d^ is the deceleration parameter. We will 
assume first that we are dealing with a canonical-scalar 
(7 > 0), so that 1 + w > (or < 0). When /i ~ 0, 
and hence H ~ const — Hini, A(cr) acts as a cosmological 
constant term, i.e. A(ct) ~ 3MpH?^f = Ainf. A major 
difference from the case of a cosmological constant Aq, 
for which pA — —PA and thus wa = —1, is that the field 
is evolving towards an analytic minimum at A((t) ~ 0. 

■ 2 ■■ 

The Gauss-Bonnet term, Q = 24^^, changes sign be- 
tween accelerating (a > 0) and decelerating (a < 0) 
solutions. In our model, the effective potential A(cr) 
(= V{a) — f{<y)Q) has to be non-negative in order that 
the vacuum becomes a possible ground state. Under the 
condition (jj)), it is required that h> —3 and consequently 
the EOS parameter w < 1. This last inequality may be 
saturated only by stiff matter for which w = 1 (and the 
velocity of sound approaches the velocity of light). In 
general, however, w <1 holds and hence A(cr) > 0. 

Under the condition (7) alone, the field potential V{a) 
may not always be positive. Then one may ask: Is there 
any condition for the (semi-)positivity of the potential? 
Below we will give the condition for V{(j) > 0, though it 
may not be required for obtaining inflationary type solu- 
tions in the present context. Furthermore the condition 
(7) fixes only one out of two arbitrary parameters in the 
model. The last arbitrary parameter may be fixed ei- 
ther by allowing one of the field variables to take a fixed 



(but arbitrary) value or by making an appropriate ansatz 
for one of the field variables. Here we consider two well 
motivated examples. 

If we use the approximation h = H/H^ = H' /H ~ 
const = ho, then, by solving |7J, we find 

^^g/Mr, ^^^g/io')^ M = u^e"-" + U2e"+'', (9) 
where Hq and Ui are integration constants, and 

a±^^(^ho~5± y^9/i2 + 54/10 + 49^ . (10) 

It is the coupling (o^-) which is of greater interest 
for ?7 ^ (77 < 0). The scale factor is given by a{t) = 
ao(l — cohot)~^/'^° , where cqHo < 0, implying that the 
universe accelerates in the range 

-1 < /z o < > g > -1. (11) 

We then wish to consider the case where h is dynamical. 
To this aim, we make no assumptions about the form of 
h, and instead consider the ansatz that during a given 
epoch we may make the approximation 

u = 8K^f{c7)H^ « ut.^ee"*""^'' (12) 

where the time can be early or late, where |uear/j/| > 
|u/ate|, and aeariy < ctiate- This ansatz IS well motivated, 
since in the context of string theory the coupling /(cr) is 
generally a sum of exponential terms (as functions of a 
or the number of e-folds) and so is the Hubble parameter. 
Interestingly, the ansatz H12|) for the form of u will allow 
us to write h in a closed form: 

H = Ho cosh l3{ri~T]^)e-^'^, (13) 
h = -(3 + ptanhPi?]- ri^), (14) 

where 77-^ is a free parameter, 

- Q; 1 J 

/3 = 4+-, /3ee -V9a2 + 72a-^ 208, (15) 

and a = Ofoariy or aiato- Writing the expressions for x(ry) 
and y{ri) is straightforward. Analogous to an inflationary 
type solution induced by a conformal- anomaly |2l| . the 
solutions given above are singularity-free. Note, for Arj = 
ri — 7]^ > 0, we always have A(cr) > for a < 1 (cf Fig.^l. 

Next, consider that u(cr) = f{a)H^ ~ const = \/a 

and X = ^^^"^-^ = Xp, where |d| > 1 and < < 1, 
which are reasonably good approximations at low energy. 
In this case, ry — ±-^7/2x~ (cr/A/p;)-fconst. The solution 
is again given by H13() - (|14|l but now 

/3 = -l±^, /3=iV(l + ")'+4xoa. (16) 
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FIG. 1: The A?; - a phase space, showing regions with V{a) > 
and V{a) < 0. The almost horizontal line with h — —3 
separates the regions between A(cr) < and A{a) > 0. Here 
we have set uume = 0.2. For a large value of uume (> 0.2), 
the regions I and IV, both with V{a) > 0, merge to a single 
region, while for a small value of uume (< 0.2), the region I 
moves to a more negative value of At;, while the region III 
moves to a more positive value of Arj (and a). In any case, the 
region with A(a) < (or /i + 3 < 0) is physically less relevant 
as it corresponds to the case where the EOS parameter w > 1. 

In particular, for xq ^ l/(4d), or /? ~ — /3, one has H oc 
Ha[l + e2/3(''-''i)]; the Hubble expansion rate decreases 
vifith 77 when Pif] — rj-^) < 0. The scalar potential and the 
coupled GB term may be given by 

V{a) = Af|,i72 ( -- - — tanh/3a-a;o ) , (17) 
\ a a J 

f{a)g = MlH^ (^±^ + A tanh/Ja) , (18) 



where a — ^J^/2x^^ n{a — cr^). Note that here a (as well 
as a) can take a negative value, so that V{(j) > for 
3 + xqch < 6/3, or, more precisely, for /i < (3 — xo)a/6. 
Initially, V{a) - f{a)g + MlH'^{i + h), but after a cer- 
tain number of e- folds, namely \P{rj — iri-^)\ > 2, so that 
|tanh/3CT| 1, the relation as above (or Eq. (O) may 
be attained by satisfying (1 + Sx^a) = (49 + 8d)(/3 — $). 
That is, the condition like (Q puts an additional con- 
straint on the model parameters but it may not exhaust 
the basic characteristics of the model. 

A question may be raised as to whether inflation is 
due to the new term in the action, f{a)Q, or sim- 
ply due to the potential V{cr). In string theory con- 
text, the coupling /(cr) may be expanded as /(cr) cx 
cosh(cr) — ln(2)] + const. At early times, fi < 0, 
the potential V{(t) is expected to dominate the Gauss- 
Bonnet contribution, f{(j)Q. At the start of inflation, 

Q = 2A (1)^ f ^ 0. Inflation is mainly due to the po- 
tential which is not essentially flat; in such a case the 
Hubble expansion rate could naturally drop to a signifi- 



cantly low value after a sufficient number of e-folds. This 
can be understood also from some earlier studies on the 
subject ^jH^I- With V{(t) = 0, the period of accelera- 
tion is short and the number of e-folds is order of unity; 
a model with vanishing potential may not be suitable for 
early inflation. At late times, the potential drops rela- 
tively faster as compared to the coupled Gauss-Bonnet 
term, so the contribution of the potential nearly equates 
or slightly exceeds the contribution of the GB term. 

Let us discuss the result (fT^ - ((T^ in some details. The 
universe starts to expand as 77 increases from rji, but it 
accelerates only when rj > rj^. Eventually when 77 > 
T]^ + 2.5/(3, the scalar field begins to freeze in, so that 
w < —1/3; the actual value of w depends on the value of 
a (= Qfeariy), scc Figs. |21 and 13 After a certain number 
of e-folds, say N, our approximation, (|12() . with the time 
being early breaks down. Sometime later, subsequent 
evolution will be controlled by (|12|) with the time being 
late. For 77 < 77/ate the universe is in a deceleration phase 
which implies that inflation must have stopped during 
the intermediate epoch. As 77 crosses rjiate, the universe 
begins to accelerate for the second time. 
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FIG. 2: The w — w' phase space occupied by the scalar 
field, in the range > 77 — 77^ > 5. From left to right a = 
-6,-5.38,-5,-4 (solid lines) and a = -hi, -0.0148, -1, -2 
(dashed lines). The universe is not accelerating for — 4 < a < 
—2. The variation of w' w.r.t. w is significant in the range 
^ ^/il ~ 2.5//3. This happens when the red-shift factor 
z < 0.8 (2.49), since /? < 4.25 (> 2), if |7?-77j is to be related 
to z via 1-1-2 = e''~''i. For |7? — ?7-^| > 2.5//3, the field is almost 
frozen, w ~ const; w < —1/3 if 77 > t;-^, while 77; > if 77 < 77-^ . 



A large shift in Hubble rate may be required to explain 
relaxation of vacuum energy from a natural value A(cr) ^ 
M^H^ ^ 10"^M| {before inflation) to a sufficiently low 
value A(cr) - M^H^ - W'^^^Ml, {after inflation), viz 
Hb ^ 10^^ eV and Ha ^ 10"^^ eV". One thus calculates 
the change in H during an inflationary epoch of N e- 
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FIG. 3: The h — H phase space occupied by the scalar 
field, in the range —5 < r] — rj^ < 5. From top to bottom 
a = -7,-6,-5,-4,-3 (solid lines) and a = 2,1,0,-1,-2 
(dashed lines). For a canonical scalar (i.e. 7 > 0), a is only 
well defined for — 6 < a < 1. The Hubble parameter decreases 
with proper time for the coupling —6 < a < 1, while it in- 
creases with proper time for the coupling q<— 6ora>l, 
heading to a value w < —1, since h> Q. In this last case, we 
get a negative kinetic energy for the field a, in some regions 
of field space, which then requires 7 < (phantom field). 



folds, by considering the ratio 



^ H{7^, + N) 



cosh/3{Ar] + N) e" 
cosh /3( At]) 



-I3N 



(19) 



where A?/ = rji ~ rj-^, with rji being the time at on-set 
of inflation. In our model, the total number of e-folds 
N required to get a small ratio, like e ~ 10~^^, depends 
on the value of a, which is related to the slope of the 
potential coming purely from the Gauss-Bonnet coupling, 
/((t) g. The value of TV would be small, N - 129, for 
a > —2 (or a < —4), while it would be large for a < 1 
(or a > -6), A^- 0(300). 

The solution l|(^- H14(l may also be used to explain a 
cosmic acceleration of the universe at late times, i.e., at 
an energy scale many orders of magnitude lower from 
inflation. The present value of the parameter h must be 
determined by observation. For jjy — rj-^l > 2.5//3, h{ri) 
varies with the logarithmic time, 77, almost negligibly. 
The equation of state parameter for the dark energy is 
w = —{2h + 3)/3. From this, and using 114() . we get 
w > —1 and (3 > (3 for —6 < a < 1, so that H decreases 
with proper time (Fig. |3Jl. For instance, if the observed 
value of the deceleration parameter corresponds to 



gobs = -1 - Khs ~ -/?-/? ~ -0.6, 



(20) 



then a — —5.3851 or —0.0148 and hence w ~ —0.73. One 
may achieve w < — 1, for 7/ 3> rj^ -I- 2.5//3, by allowing 
a < —6 or a > 1. In this case, however, since P > 
H increases with proper time. There are arguments that 



the present cosmological data may even favor a value 
w < — F22J; if this is the case, in our model, one has to 
allow a > 1 or a < —6. 

The equation of state parameter w changes signifl- 
cantly with the exponent a but not much with ry, ex- 
cept in the range \ri — r]^\ < 2.5/(3. That is, the scalar 
fleld comes to dominate the universe, presumably for the 
second time, when < 77 — ryiato < 2.5//3, and it grad- 
ually becomes a constant value, w ~ const, when ry > 
?7iatc + 2.5/(3. At late times, one finds — 1 < w < —1/3, 
depending upon the value of a, see Fig. O Whether the 
scalar field a dominates the future evolution of the uni- 
verse depends on energy density of matter system. 

In a known form for the GB coupling, derived from the 
hetcrotic string theory, ^(cr) ~ ln2 — ^cosha [isl l23|. 
It follows that /(cr) - Co + ci (e'^ +e"'^). For cr > 0, 
only the first two terms are relevant. To this end, we 
slightly modify the ansatz for u(cr), namely u(a) = uo(l+ 

e"^) with ar; < 0, which may be appropriate at some 
intermediate epoch. Note that R/2k'^ = 3{2H^ + H)/k'^ 
and f{a)g = 3u{H^ + H)/k^. We demand < |u| < 2, 
so that the GB contribution to the action is only sub- 
dominant. The solution for h is modified as 
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(3 + (3 



^1 + u^c°"i K{'q) 



(21) 

where K{ri) = {Lp^ + cLq^)/{Lp_ + cLq_) with Lp^ 
(Lq^) being associated Legendre functions of first (sec- 

, ^1 -I- u^e"''^ , and c an 
77^1 ^ and a < 0, so 
13(77 - ri^) ~ 



± 



ond) kind, Lp(Q) ^ ^ 
integration constant. For 
u{a) Mo, we find h ~ - 
—0.3944 and hence wde - 



1 2\/l3 



W 

-4 + \/l3tanh^ 
t -0.7370. This value of w is 
realized only at the asymptotic future, where u{a) at- 
tains a small (constant) value, while its present value, 
Wp)^, may be less than wde- In the presence of matter 
coupling, rim — 0.3, one defines an effective EOS param- 
eter, Weff, whose value can be greater than w^^. 

It is generally valid that the transformation from the 
string frame to the Einstein frame introduces couplings 
between standard model fields and massless degrees of 
freedom, such as, dilaton (p (a spin-0 partner of spin-2 
graviton) Simple arguments such as the absence of 
ghost, thereby guaranteeing the stability of the field the- 
ory, would suffice to rule out a wide class of scalar-tensor 
models. In our model, we have assumed, rather implic- 
itly, that dilaton is either constant or it is extremized, 
^^j^ — 0, so the effects coming from terms like X{(f>)g 
are negligible. Only the dynamical field in the model 
is the modulus cr, which is expected to have a non-zero 
mass, i.e., a non-zero vacuum expectation value. In string 
theory context 0| , there is no coupling between the field 
a and the Einstein-Hilbert term. And there is no direct 
coupling of the field a to matter, in view of precise veri- 
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fications of the weak equivalence principle. Nevertheless, 
it is useful to assume that the running of the field a is 
negligibly small at late times, namely [251 

H^a'^ < 10-^ (22) 
da 

so as to avoid a possible conflict with the time-variation 
of Newton's constant under Newtonian approximation. 

In conclusion, we find that the dark energy hypoth- 
esis fits into a low energy gravitational action where a 
scalar field is coupled to the curvature squared terms in 
the Gauss-Bonnet combination. It is established that a 
GB-scalar coupling can play an important and interesting 
role in explaining both the early and late-time accelera- 
tion of the universe, with singularity- free solutions. An 
important feature of our solutions is also that the effec- 
tive potential (or dark energy) can naturally relax to a 
sufficiently low value, A - 10"i^"Af|,, after > 125 e-folds 
of expansion. 

It is shown that one may cross the cosmological con- 
stant/phantom barrier, w — ~l, only by allowing the 
parameter a, coming from the GB coupling, such that 
the Hubble rate grows with proper time, or by allowing 
a non-canonical scalar, 7 < (phantom field). 

A number of follow-up studies may be devised includ- 
ing: a fit to determine how well the potential may be 
tuned to the present value of the cosmological constant, 
and how well slow roll parameters fit the spectral in- 
dex in an imprint on the cosmic microwave background 
anisotropy and mass power spectrum. Some of the re- 
sults appear elsewhere It is possible to generalize 
the model studied here to a system of two scalar fields, 
by introducing a dilaton field (p (281] . 
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